Towards a relativistic statistical theory 
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In special relativity the mathematical expressions, defining physical observables as the momentum, 
the energy etc, emerge as one parameter (light speed) continuous deformations of the corresponding 
ones of the classical physics. Here, we show that the special relativity imposes a proper one parameter 
continuous deformation also to the expression of the classical Boltzmann-Gibbs- Shannon entropy. 
The obtained relativistic entropy permits to construct a coherent and selfconsistent relativistic 
statistical theory [Phys. Rev. E 66, 056125 (2002); Phys. Rev. E 72, 036108 (2005)], preserving 
the main features (maximum entropy principle, thermodynamic stability, Lesche stability, continuity, 
symmetry, expansivity, decisivity, etc.) of the classical statistical theory, which is recovered in the 
classical limit. The predicted distribution function is a one-parameter continuous deformation of 
the classical Maxwell-Boltzmann distribution and has a simple analytic form, showing power law 
tails in accordance with the experimental evidence. 

PACS numbers: 05.90.+m, 05.20.-y, 03.30.+P 



I. INTRODUCTION 

In the ordinary relativistic statistical mechanics the 
entropy is assumed to have the same form of the en- 
tropy of classical statistical mechanics. Then, starting 
from the Boltzmann-Gibbs-Shannon entropy, the maxi- 
mum entropy principle yields an exponential distribution 
which exactly reproduces the Maxwell-Boltzmann distri- 
bution of classical statistical mechanics in the rest frame. 
Unfortunately the Maxwell-Boltzman distribution fails to 
explain the spectrum of cosmic rays, which undoubtedly 
represent the most famous and important relativistic par- 
ticle system 0, ■ This spectrum has a very large ex- 
tension (13 decades in energy and 33 decades in particle 
flux), exhibiting a power law asymptotic behavior. Sim- 
ilar power law asymptotic behaviors of the distribution 
function have been observed also in other relativistic sys- 
tems including plasmas (sl and multiparticle production 
processes 4]. Hence for relativistic systems, the experi- 
mental evidence suggests a non-exponential distribution 
function with power law tails. In the last 40 years, em- 
pirical non-exponential distributions with power law tails 
have been systematically used in high energy plasmas. 

A non-exponential distribution can be originated ex- 
clusively by an entropy manifestly different from the 
Boltzmann-Gibbs-Shannon one. In order to propose the- 
ories based on non-exponential statistical distributions, 
one can follow different paths, that allow to construct, 
other coherent and selfconsistent statistical theories be- 
sides the classical statistical mechanics ^SkiL^^ 

mm 

We recall that any relativistic formula, defining a phys- 
ical quantity, e.g. momentum, energy, etc, can be viewed 
as a one parameter (light speed) generalization or defor- 



mation of the corresponding classical formula. Conse- 
quently, it is natural to ask if the entropy of a relativistic 
system could be a one parameter generalization of the 
Boltzmann-Gibbs-Shannon entropy as well. 

The main goal of the present contribution is to show 
that the one-particle relativistic dynamics imposes a 
one parameter generalization of the Boltzmann-Gibbs- 
Shannon entropy. The new entropy permits the construc- 
tion of a statistical mechanics which efficiently describes 
the relativistic many-particle systems. 



II. COMPOSITION LAWS IN SPECIAL 
RELATIVITY 

We introduce the dimensionless variables momentum 
q, velocity u and total energy £ in place of the corre- 
sponding physical variables p, v and E through 



V 

u 




(1) 



being c the light speed and k a dimensionless param- 
eter. Furthermore we impose for the velocity that 
limc^oo, K^o "t"* = v'^ with v'^ = finite, in order to 
preserve the validity of the dimensionless variable def- 
initions, also in the classical limit. Within the special 
relativity it is immediate to express the deformation ef- 
fects, due to the finite value of the light speed, in terms 
of the dimensionless parameter k: 
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Consider in the one-dimensional frame S two identical 
particles of rest mass m. We suppose that the first parti- 
cle moves towards right with velocity ui while the second 
particle moves towards left with velocity 1*2. The mo- 
menta of the two particles are given hy qi = q (ui) and 
q2 = q ("2) respectively, where q (u) is defined through 
Eq. 10) ■ The energies of the two particles are given by 
£1 ^ £ {ui) and £2 = £ (^2) respectively, with £ (u) de- 
fined by Eq. We consider now the same particles in a 
new frame S ' which moves at constant speed U2 towards 
left with respect to the frame S. In this new frame, which 
is the rest frame of the second particle, we have that the 

K 

two particles have velocities given hy u'l — ui (B U2 and 



respectively, being 



Ui ® U2 



Ui + U2 
1 + k'^UiU2 



(8) 



the well known rclativistic additivity law for the dimen- 
sionless velocities. In the same frame S ' the particle mo- 
menta are given hy q'l = q{u[) and q'2 = 0, respectively. 
Analogously in S' the particle energies of the two par- 
ticle in S' results £[ = £ {u'l) and £'2 — A very 
interesting result follows from the relation £'2 = I/k^ re- 
garding the physical meaning of the deformation param- 
eter K. It is evident that the quantity I/k^ represents 
the dimensionless rest energy of the rclativistic particle. 

Let us pose now the following questions: is it possible 
to obtain the value of the momentum q ^ (of the energy 
£'i) starting directly from the values of the momenta qi 
and q2 (of the energies £ 1 and f 2) in the frame S. The 
answers to the above questions are affirmative. 

First we consider the case of the rclativistic momentum 
q'-^ and write it as q'l — q{u'i) = q{ui®U2). In ref. pTj 
it has been shown that q'l — qi ® q2 being 
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92 =qi\ll + ti'^ql + q2\l + n'^ql 



(9) 



the K-sum of rclativistic momenta. In words, the rcla- 
tivistic momentum q [ of the first particle, in the rest 
frame of the second particle S' , is the K-deformed sum of 
the momenta qi and q2 of the two particles in the frame 
S. The K-sum of the rclativistic momenta and the rcla- 
tivistic sum of the velocities are intimately related, and 
reduce both to the standard sum as the velocity c ap- 
proaches to infinity (or equivalently the parameter k ap- 
proaches to zero). The deformations in the above sums, 
in both cases, are rclativistic effects and are originated 
from the fact that c has a finite value. 

We consider now the total energy £ [ of the first par- 
ticle in the frame S'. Clearly it results £[ = ^(""'i) = 

£ (^ui © U2^ ■ In ref. it is shown that one can calcu- 
late £[ starting directly from the values, in the frame 



S, of the total energies £1 and £2- It is obtained 

£'i — £1 (B £2 where the composition law of the dimen- 
sionless relativistic total energies is defined through 

1 



£i®£2 = K''£i£2 + — J{K^£i-l){K^£^~l) . (10) 

In the following we summarize the relationships linking 
the K-sums defined through Eqs. ((HJl, © and (fTUIl 



qiui)®q{u2)=q[ui®U2j , (11) 

£iui)®£{u2)=£{ui(Su2^ , (12) 

u(gi) e u(92) = w (91 e 92) , (13) 

£{qi)®£{q2)=£{qi®q2) , (14) 

q{£i) (B q{£2) ^ q {£1 ® £2) , (15) 

u{£i) e u{£2) = u 1 e £2) ■ (16) 

Note that in the Ihs and rhs of any of the above equa- 
tions appear two different composition laws. For instance 
in Eq. Hill) in the Ihs appears the composition law of the 
dimensionless momenta ^ while in the rhs appears the 
composition law of the dimensionless velocities (jSJ. We 
emphasize that the three K-sums given by Eqs. ijSJ), © 
and H10() . emerge only when we change the frame of ob- 
servation of the particle from S to S' . 

III. LORENTZ INVARIANT INTEGRATION 

Within the special relativity a central role is played by 
the four dimension Lorentz invariant integral 



Irei =A d^p 0{po) Sip^'p^ - m'c') F 



(17) 



with A an arbitrary constant and F an arbitrary function 
depending on p = \p\. After recalling that = {E/c, p ) 
and E — ^ rn^c^ + p'^c^, one can reduce the above inte- 
gral as follows 



:F-- 



y^T+K^ Jo y^T+K^ 



dq 



■Aiiq^F , (18) 



where the constant A is fixed properly and the dimen- 
sionless momentum is introduced. We remark that in 
(|17|l the integration element d'^p is a scalar because the 
Jacobian of the Lorentz transformation is equal to unity, 
so that Irei transforms as F. Then in (|18() the integration 
element d^q/ \/l -|- K^g^ is a scalar. After introducing the 
K-differential through 



d{K)q = 



dq 



(19) 
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the integral Irei assumes the form 



Jo 



(20) 



One immediately observes that the Lorentz invariant in- 
tegral can be obtained by deforming the ordinary (clas- 
sical) integral 



Irl = 



dq Att q^ F 



(21) 



The deformation is obtained by performing the substitu- 
tion dq d^i^yq. 

It is important to note that the K-differential l|19|) can 
be obtained directly from the additivity law of the di- 
mensionless relativistic momenta © according to 



d-lnyq = (q + dq) Q q 



(22) 



IV. 



PHYSICAL MEANING OF THE 
K-DIFFERENTIAL 



We consider now two identical relativistic particles 
with rest mass m and velocities Vi and V2 respectively, 
in the frame S. The modulus of the relative velo city 
V = T^(vi,V2) of the two particles, given in ref. 
(page 20), can be written also in the form 



l^(vi, V2) 
with 



(vi e V2) 



1 



Vi X V2 
1 — ViV2/c^ 



vi e V2 



Vl - V2 
1 - ViV2/c2 



(23) 



(24) 



We perform now the calculation of the modulus of 
relative momentum of the two particles according to 
PiV) = mVl^X - V^jc^. Clearly, it results that 
Piy) = P(vi,V2) and then P(y) = P(pi,pi). At 
this point we introduce the dimensionless variables q — 
Tp/mVir and Q = Pjmv^,. After some straightforward cal- 
culation, one arrives to the following expression for the 
modulus of the dimensionless relative momentum (the 
modulus of dimensionless momentum of the particle 1 
(2) in the rest frame of other particle) 



Q(qi,q2) 



\ 



(qi e q2) - (qi X q2) 
1 + k4 (qi X Oaf 



with 



qi e q2 = qi\/ i+f^^qi - q2\/i+K2q^ 



(25) 



(26) 



The expression of the relative momentum given by Eq. 
(I25|l has been obtained in (there the formula contains 
a typing error). 



In the following, in order to explain the physical mean- 
ing of the K-differential and the K-derivative in special 
relativity, we consider that the two particles move along 
the same line so that qi x q2 = 0. 

We suppose now that the two particles have momenta 
qi = q + dq and q2 = q respectively in the frame S. To 
calculate the modulus of the momentum Q(q+dq, q) of 
one particle in the rest frame of the other particle, after 
posing q= |q|, one obtains: 



Q(q-|-dq, q) = d^^^q 



(27) 



The physical meaning of the K-differential d^^^q immedi- 
ately follows. The infinitesimal difference dq of the mo- 
menta of the particle 1 and 2 in the frame S becomes 
d^^^q if this difference is observed in the rest frame of 
one of the two particles. The meaning of the K-derivative 



A{K}q 



dq 



(28) 



readily follows. If dj dq represents the derivative in the 
frame 5, the deformed derivative d/d^i^^q represents an 
ordinary derivative in the rest frame of one of the two 
particles. 



V. THE K-EXPONENTIAL FUNCTION 

We recall that the length of any four- vector is Lorentz 
invariant. In particular for the four- momentum we have 
the dispersion relation p^p^ = rn^c^ which can be ar- 
ranged as {E/c — p){E/c+p) = m^c^ being p — |p|. This 
latter relationship after expressing the energy in terms of 



the moment according to E = -^/r 



^c^ becomes 




1 .(29) 



The latter equation, after introducing the dimensionless 
momentum q, can be also written as 



1/k 



1 



(30) 



We remark that Eq. H30(l follows directly from the rela- 
tivistic dispersion relation. Interestingly, in the classical 
limit K ^ 0, Eq. reduces to exp(— g) eyip{q) — 1 

while the dispersion relation becomes the classical one 
W = p^ /2m being W the kinetic energy. In this way we 
obtain a direct link between the dispersion relation of free 
classical particles and the ordinary exponential function. 
In the light of the above result, we reconsider Eq. H30|l 
written in the form 



exp^^j (- 1) exp^^j (g) = 1 , 



with 



1/k 



(31) 



(32) 
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We feel that the function exp^^^{q), reproducing the or- 
dinary exponential function exp(q) = exjp^^^ (g) in the 
classical limit n —^ 0, represents a one parameter rela- 
tivistic generalization of the ordinary exponential. 

Taking into account that the ordinary exponential 
results to be eigenfunction of the ordinary derivative, 
namely (d/dq) eKp{q) = exp(q), the question to deter- 
mine the eigenfunction of the K-derivative, naturally 
arises. After some simple calculation one obtains 

-T^ exp (9) = exp (g) , (33) 

SO that the K-exponential emerges as eigenfunction of the 
K-derivative d/d^K.}<l- 

The K-exponential has the following important prop- 
erty 

e'^P{K}(9i) e'^P{K}(92) = exp^^j((ji e 92) , (34) 

K 

being qi © 92 the additivity law of the dimensionless 
relativistic momenta. 

We note that Eq. (|33|l admits as unique solution the 
K-exponential and then can be viewed as a differential 
equation defining it. Analogously Eq. 134|l is a functional 
equation univocally defining the ^-exponential function. 
In conclusion we have discussed three different ways to 
introduce the K-exponential. The first way is the more 
physical one employing the relativistic dispersion rela- 
tion. The second and third ways are more mathematical 
and employ the differential equation H33() and the func- 
tional equation respectively. 

The mathematical properties of the K-exponential can 
be found in refs. ,10., ,11, 12]. One important property of 
this function is its power low asymptotic behaviour when 
its argument tends to ±00. 

VI. THE RELATIVISTIC ENTROPY 

The Boltzmann-Gibbs-Shannon entropy of a classi- 
cal particle system is defined as the mean value of the 
minus logarithm of the distribution function, namely 
S{f) = -(ln(/)) = -Jd^pf ln(/). In order to define 
the entropy of a relativistic particle system, we employ 
the same definition with the only difference that the or- 
dinary logarithm now is replaced by the K-logarithm, de- 
fined as the inverse function of the ^-exponential, namely 



The K-logarithm is a one parameter continuous deforma- 
tion of the ordinary logarithm which is recovered in the 
limit K ^ 0. The reason of the modification of the classi- 
cal entropy is due to the fact that the ^-exponential and 
then the K-logarithm naturally emerge within the special 
relativity at the place of the ordinary exponential and 
logarithm appearing in classical statistical mechanics. 
The new entropy is given by 

Sif) = -{ In^Jf)) =-Jd'pf In^Jf) , (36) 

and after maximization under the proper relativistic con- 
strains, yields the distribution 

^ = "^"P<m(^ Xk^f ) ' (37) 

with a = [(1 — k)/(1 -I- k)]^^'^'^ and A = ^/l — k^. Being 
both these parameters real, |«;| < 1 follows. The above 
distribution results to be quite different with respect to 
the relativistic Maxwell-Boltzmann distribution where 
the K-exponential is replaced by the ordinary exponen- 
tial. The distribution (|37|l . in the global rest frame where 
the hydrodynamic four vector velocity is C/" = (c, 0, 0, 0) 
and in absence of external forces (A" = 0), simplifies to 

/-ex p,^>(-^) , (38) 

being W the relativistic kinetic energy. The latter dis- 
tribution presents, for W — > 00, power law tails, namely 
/ oc W~^l'^^ in contrast to the Maxwell-Gibbs distri- 
bution which decays exponentially. Moreover the distri- 
bution H38|l fits very well the experimental data of the 
cosmic ray spectrum |_1 1| . 

It is already known in the literature that the en- 
tropy (|36|l has the following properties: i) it is non- 
negative and achieves its maximum value at equiprobabil- 
ity, /(p) = 1/ri for Vp; and this value is 5 = In^^jO; ii) it 
is concave so that the system is stable in thermodynamic 
equilibrium; iii) it satisfies the Lesche stability condition 
and then it is physically meaningful; iv) it generates a 
selfconsistent and coherent statistical mechanics. Fur- 
thermore it results that: v) this statistical mechanics can 
be obtained as stationary case of a kinetics governed by a 
generalized relativistic Boltzmann equation; vi) this gen- 
eralized Boltzmann kinetics obeys the H-thcorem which 
represents the second law of thermodynamics. 
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